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SIDELOBE  LEVEL  OF  AN  ADAPTIVE 
ARRAY  USING  THE  SMI  ALGORITHM 


I.  INTRODUCTION 

An  adaptive  antenna  array  adjusts  its  antenna  element  settings  so  as  to  null  out  interfering 
sources  while  maintaining  a  beam  in  a  desired  signal  direction  [1].  Two  such  adaptive  array  confi¬ 
gurations  are  the  “full  up”  adaptive  array  and  the  sidelobe  canceller  (SLC)  illustrated  in  Figs.  1(a) 
and  1(b),  respectively.  The  full  up  adaptive  array  adjusts  all  the  weights  on  identical  antenna  array 
elements  while  only  the  auxiliary  antenna  weights  are  adjusted  on  the  SLC  configuration.  The  full  up 
adaptive  array  may  have  some  mainbeam  constraints  associated  with  it  so  that  a  desired  mainbeam  can 
be  maintained  [1]. 


PHASED  ARRAY 


Fig.  la  —  Full  up  adaptive  array 
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The  weights  on  the  full  up  adaptive  array  are  adjusted  so  as  to  maximize  the  output  signal-to- 
noise  power  ratio,  and  the  weights  on  the  SLC  are  adjusted  so  as  to  minimize  the  output  noise  power 
residue,  if  the  noise  environment  is  not  known  a  priori,  these  weight  settings  cannot  be  set  a  priori 
but  n.ust  be  estimated  from  a  finite  set  of  incoming  data  on  the  input  channels.  Thus  the  weight  set¬ 
tings  will  have  perturbations  about  the  quiescent  optimum  weight  settings.  These  perturbations  result 
in  a  rise  in  the  adaptive  array  sidelobe  level  above  the  quiescent  adaptive  array  sidelobe  level.  We 
call  the  nonquiescent  sidelobe  level  that  is  due  to  finite  sampling  and  averaging,  the  transient  sidelobe 
level. 


This  report  presents  an  analytical  result  for  predicting  this  rise  in  the  adaptive  array  sidelobe 
level  by  use  of  the  Sampled  Matrix  Inversion  (SMI).  The  SMI  algorithm  [2]  is  an  open-loop,  rapidly 
converging,  adaptive  array  implementation  whose  convergence  rate  is  independent  of  the  external 
noise  environment.  For  many  years,  it  has  been  considered  a  baseline  for  fast  converging  adaptive 
array  algorithms.  Brennan  [3]  has  presented  some  theoretical  and  computer  simulation  results  for  the 
effects  of  the  SMI  algorithm  on  the  sidelobe  level.  His  theoretical  results  pertain  to  the  single  auxili¬ 
ary  SLC.  We  generalise  these  results  by  examining  an  array  with  any  number  of  inputs.  Moreover, 
many  of  his  observations  made  as  a  result  of  computer  simulations  are  given  a  theoretical  basis. 

The  SMI  algorithm  is  briefly  reviewed  in  Section  II,  and  pertinent  SMI  theorems  are  presented 
in  Section  EH.  An  exact  expression  for  the  transient  sidelobe  level  is  developed  in  Sections  IV  and  V. 
A  discussion  of  this  result  appears  in  Section  VI,  techniques  for  lowering  the  sidelobe  level  are 
presented  in  Sections  VII  and  VIII,  and  other  techniques  that  can  be  used  to  reduce  the  transient 
sidelobe  level  are  discussed  in  Section  IX. 

D.  THE  SMI  ALGORITHM 

In  this  section,  we  briefly  review  the  SMI  algorithm  and  establish  the  notation  and  assumptions 
relevant  to  the  succeeding  development.  Applebaum  (4]  has  shown  that  the  optimal  (maximization  of 
the  average  output  signal-to-ncise  ratio)  adaptive  weighting  of  an  N  element  array  is  given  by 

w  =  nM~xs*  (1) 


where  w  is  the  optimal  weighting  vector  of  length  N,  s  is  a  normalized  steering  vector  related  to  the 
direction  of  arrival  of  the  desired  signal,  M  is  the  steady  state  covariance  matrix  of  the  inputs,  n  is  an 
arbitrary  constant,  and  *  denotes  the  complex  conjugate.  More  formally  we  write 


W  =  (w0,  wx,  ,  Wft  _  ])^" , 

(2) 

X 

II 

X 

o 

* 

1 

(3) 

and 

M  =  E\X*XT), 

(4) 

where  T  denotes  the  vector  (or  matrix)  transpose  operation,  £|) 

denotes  the  expected  value,  and  X0, 

Xx,  •  • ' ,  XN-\  are  the  inputs  of  the  adaptive  array.  Because  M  is  generally  not  known  a  priori,  it  and 
thus  the  adaptive  weights  must  be  estimated  from  sampled  data.  Hence,  time  samples  of  the  nth 
channel,  Xn(k),  are  taken  where  k  indexes  the  sampled  data.  Note  that  the  data  on  each  channel  for  a 
given  k  are  time  coincident. 
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r^r  this  development,  we  make  the  following  assumptions  unless  otherwise  noted  (the  same 
assumptions  were  made  in  Ref.  2). 

1.  The  Xlt  X2, . . . ,  XN  are  zero  mean  stationary  Gaussian  random  variables  (r.v.). 

2.  Xfl,(*i)  ‘s  independent  of  XBj(k2)  for  nt  *  n2  or  *  k2. 

3.  The  estimated  adaptive  weights  are  computed  from  an  input  data  set  that  is  independent  of 
the  data  that  are  adaptively  weighted  by  the  adaptive  array.  (We  call  this  nonconcurrent 
processing.) 

4.  The  desired  signals  are  not  present  duiing  the  adaptive  weight  computation. 

With  respect  to  Assumption  3,  if  the  weights  are  computed  and  applied  to  the  same  set  of  data,  we 
call  this  concurrent  processing  and  the  output  residue  is  called  the  concurrent  output. 

We  write  the  sampled  input  vector  as 

X(*)  -  (X0(Jt),  *,(*) . X*_,(*))r.  (5) 


For  the  SMI  algorithm,  the  covariance  matrix  A#  is  estimated  as 

(6) 

Kk"l 

where  the  caret  over  the  M  denotes  an  estimate,  K  is  the  total  number  of  independent  samples  per 
input  channel,  and  K  Ss  N  so  that  M  is  not  necessarily  a  singular  matrix.  Note  that  M  is  calculated 
in  “batch”  style  from  a  block  of  K  by  N  input  data.  Thus  after  M  is  calculated,  we  find  the  esti¬ 
mate  of  w  as 

*  =  nM~'s*.  (7) 

As  K  —  oo,  M  —  M,  so  that  for  an  infinite  number  of  samples,  the  optimal  array  weighting  vector 
is  obtained. 

In  many  instances  in  the  following  discussion  we  refer  to  the  SLC  configuration  of  an  adaptive 
array.  For  this  configuration 


s  =  (1,  0,  . . . ,  0)r  and  w0  =  1 ;  (8) 

i.e.,  the  desired  signal  is  present  only  in  the  Oth  channel,  which  we  call  the  main  channel.  The  N  —  1 
other  channels  referred  to  as  the  auxiliary  channels  are  used  to  cancel  unwanted  signals  (noise)  from 
the  main  channel.  Note  that  the  condition  that  w0  =  1  is  a  constraint  on  the  optimal  weighting  vec¬ 
tor. 


—  (Xj,  X2,  . . .  ,  Xn~i)t  , 

(9) 

w0  =  (w,,  w2, ... ,  wN_  y , 

(10) 

X/fc)  =  (*,(*),  X2(*) . Xv_,( k))T, 

(ID 
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Ma  =  £{XJXj) , 


ram=E[X*X0)t  (13) 

where  the  subscript  a  refers  to  the  auxiliary  channels.  It  can  be  shown  that  the  optimal  weighting 
vector  of  the  auxiliary  channels  for  the  SLC  configuration  is  given  by 

=  4  AC'fam  •  (14) 

Again  with  no  a  priori  information,  it  is  necessary  to  estimate  Ma  and  ram .  These  estimates  are  given 
by  the  expressions 

Ma  =  X*(*)Xj(*)  (15) 

**  =  l 


=  4  £  x*(k)x0(k). 


For  the  steady  state  weights  or  optimal  weights,  we  can  write 

N- 1 

x0  =  E  w*x*  +  e ' 

n  =  1 


where 


E[eX*a)  =0. 


Equation  (18)  is  merely  an  exemplification  of  the  orthogonality  principle. 

ffr  SMI  THEOREMS 


In  this  section,  we  present  and  prove  (or  cite  references  for)  a  variety  of  theorems  related  to  the 
SMI  algorithm.  We  use  these  theorems  in  later  sections  to  derive  closed  form  solutions  for  the  tran¬ 
sient  sidelobe  level  of  an  adaptive  array. 

Theorem  1:  If  the  input  vector  X  is  transformed  by  a  nonsingular  matrix  A ,  then  the  estimated 
transformed  weights  w'  are  related  to  the  estimated  untransformed  weights  (V  by  the  relationship 


1 


$ 


!V.» 


% 


m 


Proof:  We  know  that 


w  =  A  rw' . 


where 


w'  =  M'-'s'*  , 


M'  =  l AX(k)]*  [AX{k)]T 

K  k  =  i 
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and 


Using  Eq.  (22),  we  can  show  that 

! 


s'  «  4s. 

AT  -  A*MAt . 


(22) 

(23) 


Using  Eqs.  (21)  and  (22)  in  Eq.  (20),  we  find  that 

*'  -  l<i4*M4rj-1(j4s)*  (24) 

=  (4r)_1  M~ln* 

=  mV1*. 


The  theorem  follows  from  Eq.  (24). 

Similarly,  we  can  show 

Theorem  2:  If  the  adaptive  army  is  in  the  SLC  configuration  and  the  auxiliary  input  vector  X<,  is 
multiplied  by  a  nonsingular  matrix  transform  A ,  then  the  estimated  auxiliary  weights  *'  are  related  to 
the  estimated  untt^is  formed  auxiliary  weights  wa  by  the  relationship 

*a  =  ATfy.  (25) 

Theorem  3  is  cited  in  Ref.  2. 

Theorem  3:  There  exists  a  unitary  matrix  U  that  transforms  the  steering  vector  s  into  the  vector 
(1,0,...,  0)r,  where  s's  =  1. 

Theorem  3  is  important  to  our  development  because  it  allows  us  to  transform  any  unconstrained  adap¬ 
tive  array  into  the  SLC  configuration. 

Theorem  4:  If  the  input  data  vector  X  or  Xa  is  multiplied  by  a  unitary  matrix  U ,  then 

£{*'*}  =  £{*"*'}  (26) 

or  for  the  SLC  configuration 

=  £{*„'*').  (27) 

Equations  (26)  and  (27)  follow  directly  from  the  fact  that  UU'  =  I . 

N  - 1 

Theorem  5:  IfX 0  =  £  cnXn  and  the  adaptive  array  is  in  the  SLC  configuration,  then  the  concurrent 

n  =  1 

output  of  the  SMI  canceller  is  exactly  zero  independent  of  K ,  the  number  of  independent  samples  per 
channel  taken. . 


5 


KARL  GERLACH 


Proof:  We  write 

*o(*>  -  "t'cMk)  -  X/(*) c,  (28) 

*  “1 

where  we  define 

c  *  (cj,  c2, ...  ,cw_l)r.  (29) 

Thus 

-  “£xS(*)X.r(*)c  (30) 

A*«i 


Hence  if  we  solve  ASf0^0  =  =  #ac,  we  see  that  *  c,  and  the  theorem  follows. 

Theorem  6  results  from  Theorem  5. 

N-l 

Theorem  6:  If  X0  —  £  cHXn  +  e,  then  the  concurrent  output  noise  residue  is  independent  of 
N- 1 

r<vv 

n  *  1 

n  (  Theorem  6  is  exemplified  by  Fig.  7.  Here,  we  implement  two  parallel  cancellers;  one  cancels 
£  c„  X„ ,  and  the  other  cancels  e  .  Note  that  separate  estimated  weighting  vectors  are  computed  for 

n  =  l 

each  canceller.  As  a  result  of  Theorem  5.  the  left-hand  canceller  always  has  zero  residue.  Hence  the 
final  lesidue  R  is  dependent  on  only  the  right-hand  canceller’s  output  residue. 


Fig.  2  —  Equivalent  SMI  canceller 
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Theorem  7  was  first  proved  by  Reed  et  ai.  [2]. 

Theorem  7:  If  assumptions  1  through  4  are  satisfied  and 

_S|W2_ 

(S/N)m  ' 

where  (S  /N)  is  the  sampled  signal  to  noise  ratio  given  by  the  expression 

($/N)  - 

and  S/N)  is  the  optimal  signal-to-noise  ratio  given  by 

(S/N)  =  w'A/w , 


then  the  probability  density  Junction  (p.d.fi)  of  z  is 

K\ 


P(z)  - 


(N  -  2)!  IK  +  1  -  ,V)! 


(1  -  z)N  ,  0  s  z  s  1. 


(3!) 


(32) 

i 

i 


(33) 


(34) 


From  this  theorem,  Brennan  and  Reed  [5)  showed  that 

Theorem  8:  If  assumptions  1  through  4  are  satisfied  and  the  adaptive  array  is  in  the  SLC  configura¬ 
tion,  then 


°SMl(^>W) 


"min 


K-N  +  \ 


K  a  N 


(35) 


where  ,  N)  is  the  average  output  noise  power  residue  of  the  N  input  SMI  canceller  using  K 

samples,  and  a^n  is  the  minimum  output  noise  power  residue  (K  =  ao). 


We  can  show  that 


and 


°s2mi(*.  N)  =  £(  |  A'0|  2J  -  £{*'Af0  *„) 

4in  “  fl  1^0 12)  ~  <MUWU. 


(36) 

(37) 


IV.  SIDELOBE  LEVEL  MEASURE 

In  this  section,  we  derive  an  expression  for  the  sidelobe  level  contribution  of  the  adaptive  array 
that  comes  as  a  result  of  estimating  the  optimal  adaptive  weights.  We  consider  the  two  adaptive  array 
configurations,  full  up  and  SLC. 

We  simplify  the  full  up  adaptive  array  analysis  by  transforming  the  full  up  array  into  an  SuC 
configuration  by  invoking  Theorem  3.  (The  appendix  gives  a  simple  procedure  for  finding  U  when 
'Jn  s  =  (1,  e i*t  e2J*,  •  •  ■ ,  e(N~l)i*)T).  Let  U  be  the  unitary  transform  such  that 
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We  place  the  SLC  constraint  (w0  "  1)  on  the  adaptive  weights  that  are  computed  after  this  transfor¬ 
mation.  We  shall  see  that  this  unitary  transformation  does  not  change  the  transient  sidelobe  level. 

Let 

gq($)  be  the  quiescent  antenna  pattern  with  no  adaptation  (0  s  4  <  2»), 
g^  (4)  be  die  quiescent  adaptive  antenna  pattern  ( K  -  »), 
g,(4>)  be  the  adaptive  antenna  pattern  (finite  K), 

be  the  quiescent  adaptive  weighting  vector  {K  ■  »,  SLC  configuration), 
w  be  the  estimated  adaptive  weighting  vector  (finite  K)t  a  random  variable, 

Aw  «  w  -  wf ,  a  random  variable, 

HO)  =  (1,  eJ ♦,  e **♦, . .. ,  eiN~X)i4)T,  and 

Um  be  (he  nth  row  of  {/. 

The  auxiliary  inputs  have  the  following  antenna  patterns: 

*„(*)  -  u/v(*),  n  -  1,  2 . N  -  1.  (39) 

We  define  a  vector  of  length  N  —  1  of  these  auxiliary  antenna  patterns: 

8  (*)  =  Of  i(*).  •  .  gN-imT.  (40) 

By  using  the  above  defined  quantities,  it  is  simple  to  show  that 

£<,(*>  =*,(*)  -  wrg(4>)  (4D 

=  gq(<t>)  -  (W„  +  Aw )TU\(<t>). 

Taking  the  expected  value  over  the  Aw  of  the  magnitude  squared  of  Eq.  (41)  results  in 

EawI  l&i(fl!>)|2)  =  \gq(<t>)  -  W Ji/V($)|2  +  E±w\  Aw't/*v*(0)vr(<^)t/rAwl.  (42) 

The  fact  that  £Aw{Aw|  =  0  was  used  in  the  above  derivation.  Note  that 

gqatt)  =  8q(<t>)  ~  (<£)  .  (43) 

so  that 


E*wi\ga(<t>)\2)  =  \gqa&)\2  +  £A,,lAw't/*v*(*)vr(*)f/rAw}. 


(44) 
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Thus  we  Me  tint  the  adaptive  transient  sidelobe  level  contribution,  A SL* ,  is  given  by 

ASX*  ~  £AwfAw,t^v*(d)vr(d)£/rAw).  (45) 

The  average  contribution  is  found  by  averaging  ASL,  over  0  s  4  <  2*.  Because 

£*l**(d)vr(d)l  -  ±  £  v*(4)vr(4)*4  (4ft) 

-  /. 

where  /  is  the  (AT  -  1)  x  (N  -  1)  identity  m-  rix,  it  follows  that 

ASZ*  -  EIAw'Aw),  (47) 

where  ASX*  is  the  average  transient  sidelobe  level  for  a  full  up  adaptive  array,  and  we  have  dropped 
the  Aw  subscript  from  the  expected  value  notation. 

Consider  the  SLC  configuration  as  shown  in  Fig.  1(b).  Here 

gM(d)-GeKV  n-1,2 . N  -  i ,  (48) 

where  mnt{0,  1 ,  ....  AT  -  1),  the  m„  are  distinct  for  each  n ,  and  G  is  an  arbitrary  gain  associated 
with  the  auxiliary  antenna  elements.  Similar  to  the  p.  reeding  development  for  the  foil  up  array,  we 
can  show  that 


ASL,,  =  G2£|Aw/Aw). 


(49) 


Hence  the  SLC  configuration's  average  sidelobe  level  expression  is  the  same  as  the  full  up  adaptive 
array  except  for  the  factor  G2.  We  use  Eq.  (49)  for  both  configurations  by  defining  G  =  1  for  the 
foil  up  array. 

Note  that_as  a  result  of  Theorem  4,  the  unitary  matrix  transformation  of  the  input  channels  does 
not  change  A SL„ .  Also  due  to  the  form  of  A SLa  given  by  Eq.  (49),  Assumption  3  as  given  in  Sec¬ 
tion  II,  need  not  hold. 

V.  SIDELOBE  LEVEL  DERIVATION 

In  this  section,  we  derive  an  exact  expression  for  the  average  transient  sidelobe  level,  A5£y,.  In 
the  derivation  we  consider  only  the  SLC  configuration,  but  the  result  is  applicable  for  the  foil  up 
array. 

Consider  the  concurrent  output  residue  of  the  SMI  canceller  as  shown  in  Fig.  3.  Here 

R  =  X0  -  *rX„  (50) 

=  Xq  -  (W,  +  Aw)rXa  . 
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AUX  CHANNELS 
*o(  \\  X, 


|  x’|  x°|  .  .VI 


Fig.  3  —  Representation  of  SMI  canceller 


Using  Eq.  (17),  we  can  write 


R  =  e  -  AwrXa  ,  (51) 

where  e  is  statistically  orthogonal  to  Xa .  In  fact,  invoking  Theorem  6,  we  see  that 

Aw  =  ,  (52) 

where 

=  T  tX*a(k)e(k)  (53) 

*  k  =  1 

and  e(k),  k  —  1 , 2,  . . . ,  K  are  the  time  samples  of  e . 

Since  the  auxiliary  input  covariance  matrix  is  hermitian,  we  can  write 

Ma  =  &  A4>,  (54) 

where  A  is  a  diagonal  matrix  of  the  eigenvalues  of  M  and  $  is  a  matrix  of  eigenvectors  of  M .  Note 
that  we  assume  that  all  powers  are  referenced  to  the  internal  noise  level,  which  we  set  equal  to  one. 
The  eigenvalues  can  then  be  divided  into  two  classes:  significant  eigenvalues  >>  1,  and  noise  eigen¬ 
values  that  are  approximately  equal  to  1  or  the  internal  noise  power  level. 

We  transform  the  auxiliary  input  by  the  matrix  transform  A~l/2  4>*  as  shown  in  Fig.  4,  where 
A~1/2  designates  a  diagonal  matrix  whose  diagonal  elements  are  equal  to  the  square  root  of  the 
corresponding  diagonal  elements  of  the  inverse  of  A.  Using  Theorem  1,  we  see  that 

Av/  =  (A~1/2$*)'r  Aw'  (55) 

=  4>'A~ 1/2  Aw'. 


Thu’ 


£{Aw'Awj  =  £{Aw''A  1/2$$'A  1/2Aw'J. 


(56) 
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Now  let 


AUX  CHANNELS 

1  V 

|  Xl  I  X2  |  XN- 1 


a‘1/2* 


SMI 

CANCELLER 


Fig.  4  —  Orthogonalization  of  the 
auxiliary  channels 


Aw'  =  (Ah»{,  A w^,  ,  Aw^_,)r , 


*o  =  (1,  Awj,  A  w$,  ....  A  w'N_x)T , 


and  Xlt  X2,  .  be  the  eigenvalues  of  Ma  (and  the  diagonal  elements  of  A).  Because  =  /, 

we  can  rewrite  Eq.  (56)  as 

E  [ Aw1  Aw)  =  £  -1-E  [  |  A<  | 2) .  (59) 

n  =  l 

Because  of  symmetry,  we  know  that 

£{|Aw{|2}  =  £{|A>^|2}  =  ...  =  £{  |  Atv^_!  | 2) .  (60) 

Also,  the  nonconcurrent  output  noise  power  residue  is  given  by 

<7smi(*.  AO  =  £{*oM'w0}>  (61) 


where  M'  is  the  covariance  matrix  of  the  inputs  (main  and  auxiliaries)  of  the  SMI  canceller  shown  in 
Fig.  3.  Because  all  the  inputs  are  independent  random  variables,  M'  is  an  N  x  N  diagonal  matrix 
with  the  first  diagonal  element  equal  to  a^in,  and  all  other  diagonal  elements  equal  to  1  so  that 

«&!<*,  AO  =  <4ta  +  E^l|A<|2}. 


(62) 
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From  Theorem  8  and  Eq.  (62)  it  follows  that 

n- l  N  —  l 

|2»  =  "  1 


K  &  N . 


From  Eqs.  (60)  and  (63),  it  follows  that 


E{  |Aw;|2J  = 


-  JV  +  1 


^min*  n  ~  1,2,  ...  ,N  1. 


Hence,  substituting  Eq.  (64)  into  Eq.  (59),  we  have  shown  that 


i.  ,  umin  1 

£(iV4w!  *  K~N  +  T  x7 


or  the  average  transient  sidelobe  level  contribution  is  given  by 

G20l  N-  1  1 

w  vnun  v-.  1 


A5La  =  ^ 


-  N  +  1  ",  X„  ’ 


K  ±N. 


(63) 


(64) 


(65) 


(66) 


VI.  DISCUSSION 

We  see  from  the  expression  derived  for  the  average  transient  sidelobe  level,  Eq.  (66),  that_ASLa 
is  eigenvalue  dependent:  the  more  significant  eigenvalues  (X>>  1)  there  are,  the  smaller  A SLa  is. 
In  addition,  we  observe  that  ASLa  increases  with  the  gain  on  the  auxiliary  elements  G .  As  Brennan 
[3],  pointed  out,  there  is  a  tradeoff  in  selecting  G  in  that  although  A SLa  increases  with  G,  the  output 
noise  residue  after  cancellation  decreases  because  the  adaptive  weights  are  smal'er,  resulting  in  less 
amplification  of  the  internal  noise  in  the  auxiliary  channels.  The  effect  of  G2  on  the  eigenvalues  is  to 
increase  the  significant  eigenvalues  proportionally  (those  » 1)  but  leave  the  noise  eigenvalues 
unchanged  and  approximately  equal  to  one.  Finally,  we  see  from  Eq.  (56)  that  as  expected, 
ASLa  —  0  as  K  —  oo ,  and  that  the  transient  sidelobe  level  increases  as  the  steady  state  residue  a2in 
increases. 


We  note  that  for  the  full  up  adaptive  array  the  eigenvalues  we  are  referring  to  are  those  of  the 
N  —  1  auxiliary  inputs.  Hence  these  depend  on  the  unitary  matrix  transform  U ,  which  created  those 
auxiliary  inputs.  However,  since  the  number  of  significant  eigenvalues  is  generally  equal  to  the 
number  of  degrees  of  freedom  (DOFs),  Woof,  needed  to  effectively  suppress  a  given  external 
environment  noise,  we  see  that  the  number  of  significant  eigenvalues  is  a  constant.  Hence,  if 
Xj  >  X2  >  . . .  X/vTOp  »  1  and  X„  ~  1  for  N  >  Ndq f,  then  we  can  approximate  ASLa  as 


A  SLa 


r  2„2 

G  Omi 


N  ~  Ndqf  —  ^ 
K  -  N  +  1 


(67) 


Equation  (67)  can  be  used  to  properly  specify  the  number  of  independent  samples  per  channel  K 
necessary  for  the  adaptive  array  sidelobes  to  settle  within  some  arbitrary_value  of  the  quiescent  adap¬ 
tive  sidelobe  level.  If  SLa  is  the  average  sidelobe  level  of  ga((f>)  and  SLqa  is  the  average  quiescent 
sidelobe  level  of,  gqa  (<)>),  then  we  can  write 
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SLa 


SLqa 


+  G24in 


N  —  N  pop  ~  1 
K  -  N~+  1 


(68) 


We  find  the  AT  to  be  such  that  SLa  is  within  3  dB  of  the  quiescent  adaptive  sidelobe  level,  SL^ .  This 
is  found  simply  to  be 


*3dB  -  (N  -  1) 


G20  min 
nun 


1  +  — =^r 


SL. 


'qa 


G2oL 


~  ^DOF  — 


SL. 


qa 


(69) 


If  SL is  the  average  quiescent  sidelobe  level  when  the  sidelobe  cancellor  is  disabled  (all  weights 
set  equal  to  zero)i_then  normally  the  adaptive  array  is  designedjo  that  the  quiescent  adaptive  sidelobe 
level  is  equal  to  SLq  and  =1.  If  we  define  Gm  =  G2/SLq  to  be  the  gain  margin  of  the  auxili¬ 
ary  antenna  elements  with  respect  to  the  main  antenna’s  average  sidelobe,  then  we  can  show  by  using 
Eq.  (69)  that 


K2dB  =  (N  -  1)  (1  +  Gm)  -  MdqfGa,-  • 


(70) 


When  there  is  no  external  interference,  then  SLa  is  at  a  maximum.  For  this  worst  case, 
A! dof  =  0,  so  that 


SLa  =s  G2 


N  -  1 

k  -  n  +  r 


(71) 


In  addition,  for  this  worst  case,  KMS  is  at  a  maximum,  SLq  is  exactly  equal  to  SLqa ,  and  a^in  =  i. 
We  use  this  worst  case  scenario  to  specify  K3dB: 


K3 dB  -  W  -  1)  (1  +  Gm). 


(72) 


We  set  Naux  —  N  -  1,  the  number  of  auxiliary  input  channels,  and  plot  KMB/Naux  vs  Gm  in  Fig.  5. 
Note  that  for  Gm  » 1,  A3dB  is  directly  proportional  to  the  gain  margin.  Hence  large  gain  margins 
result  in  long  settling  times  for  the  adaptive  sidelobes. 

VII,  SIDELOBE  LEVEL  REDUCTION  BY  NOISE  INJECTION 

In  this  section,  we  give  the  theoretical  basis  to  a  technique  discussed  in  Refs.  3,  6,  and  7,  which 
significantly  reduces  the  settling  times  of  the  adaptive  sidelobes.  For  this  technique,  arbitrary 
independent  noise  is  injected  into  each  auxiliary  channel  (SLC  configuration).  This  noise  is  indepen¬ 
dent  from  sample-to-sample,  channel-to-channel,  and  has  noise  power  equal  to  a2.  It  can  be  shown 
that  in  the  steady  state  the  adaptive  auxiliary  weights  are  given  by 

=  (Ma  +  cr2/) —  1rma  .  (73) 

After  the  auxiliary  weights  have  been  calculated  by  using  the  SMI  algorithm,  these  weights  are  then 
applied  to  auxiliary  data  (concurrent  or  nonconcurrent)  that  do  not  contain  the  injected  noise.  The 
noise  power  of  the  injected  noise,  a2,  is  chosen  so  that  the  steady  state  noise  power  residue  is  not 
significantly  increased  (normally  within  1  dB).  Brennan  [3],  using  computer  simulation,  showed  that 
o2  can  be  chosen  to  be  much  greater  than  one  (the  referenced  internal  noise  power  level)  without 
seriously  degrading  cancellation  performance. 
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MAIN  AUXILIARIES 


Fig.  6  —  SMI  canceller  using  noise  injection 


Again  if  we  consider  the  worst  case  (no  external  noise),  then 


^3dB 

N  ~  _2  ' 

l'aux  O 


Vm.  SETTLING  TIME  REDUCTION  BY  NOISE  INJECTION 


(77) 


If  we  use  the  noise  injection  technique  described  in  the  preceding  section,  we  can  also  decrease 
the  seti)in<>  time  of  the  nonconcurrent  output  noise  power  residue.  To  see  this  we  define  the  follow¬ 
ing  quantities: 


r  lw(5.5i).  E{  \R  |2]  T{  is  the  steady  state  (K  =  oe)  output  noise  power  residue  by  use  of  the 
noise  injection  technique,  T) 

E{  |  R  | 2)  T{  is  the  output  noise  power  residue  for  finite  sampling  by  use  of  the  noise  injection  tech¬ 
nique,  T) 

E{  |  Rin  |2)  T{  is  the  steady  state  output  noise  power  residue  if  weights  are  applied  to  auxiliaries  with 
injected  noise,  T) 

E{  |  Rin  | 2)  T(  is  the  output  noise  power  residue  for  finite  sampling  is  the  weights  are  applied  to  auxi¬ 
liaries  with  injected  noise,  T) 

wa,in  T(  is  the  steady  state  auxiliary  weights  by  use  of  the  noise  injection  technique,  and  T) 

$/a  T{  is  the  estimated  auxiliary  weights  by  use  of  the  noise  ii^ection  technique.  T) 


Aw  a, it  ^a,in  '  w a, in 


(78) 


A fa, in  =  K  + 


(79) 
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Now  we  know  that 


E{\R\2}  w-rttoXa|a), 


and  we  can  show  that 


EIISI1!  -  £t|*o-*.r*S,  Is! 


=  £{|^0  -  (wa +  Aw aJ„y\a  |2} 


=  £(|/?  |2}  +  E { Aw '  Ma  Awa  } 


<i 

i,fJ 

|| 

R. 

Fv 


=  E[\R  I2}  +  £tAw'inA#a>/llAwa<IJ  +  o2E{ Aw'^Aw^J 


Using  Theorem  8,  we  can  show  that 


f  ~  1 

*(l*m  I2)  +  |2) 


£  l ,/«M«  Awa  .,„)=£{  |  Km  | 2j 


2)  AL-1 


AT-  N  +  1 


In  Section  V,  we  show  that 


r>  i  a  I  a  >  I  )  1 

£{Aw'>inAwa>jJ  =  - - -  ^  . 

K  -  N  +  l  Hml  \n  +  a2 

Hence,  substituting  Eqs.  (83)  and  (84)  into  Eq.  (81)  and  normalizing. 


E<  I 

I‘t| 

Pi 

II 

II 


Now  normally  a2  is  chosen  so  that  E\  ( Rin  |  2]  —  £(  I R  1 2).  We  then  rewrite  Eq.  (85)  as 


c  l  |a  |  }  K-N+l  K  -  N  +  1  X„  +  o* 

We  see  from  Eq.  (86)  that  the  noise  eigenvalues  (X  =  1)  in  the  summation  will  dominate,  and  if 
a2  » 1 ,  then 

I2)  K  N  -  \  -  N dof 


E[\R  |21  K  -  N  +  1  K  -  N  +  1  ’ 


1 

I 

I  tl 


E[\R  l2}  K  -  N  +  N pop  4-  1 
E\\R  |21  K  -  N  +  1 
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If  Vi  calculate  how  many  samples  are  needed  so  that  the  output  noise  power  residue  is  within  3 
d£  of  the  quiescent  value,  we  find 


^3dB  =  N  —  1  +  A^DOF* 


(89) 


Note  that,  N  s  K 3^  <  IN ,  depending  on  the  number  of  DOFs  used.  Also,  Eq.  (89)  implies  that 
the  less  DOFs  used,  the  faster  the  settling  time.  With  no  noise  injection,  =  2N  -2  [3]. 
Hence  noise  injection  decreases  the  adaptive  algorithm  settling  time. 

Another  method  to  reduce  the  settling  time  of  the  SMI  algorithm  is  to  estimate  the  number  of 
DOFs  needed  and  set  the  number  of  auxiliaries  equal  to  N^qf-  One  method  for  estimating  Nxxyr  *s  to 
count  the  number  of  significant  eigenvalues  of  the  estimated  auxiliary  covariance  matrix.  A  less  com¬ 
putationally  intensive  technique  is  to  use  the  open-loop  Gram-Schmidt  (GS)  canceller  [8,9]  to  imple¬ 
ment  the  SMI.  Here,  the  noise  power  at  each  level  of  the  main  channel  in  the  GS  structure  is  moni¬ 
tored  by  using  finite  averaging.  Note  that  for  a  finite  K  that  if  N  >  Adof>  the  output  noise  power  in 
the  main  channel  decreases  and  then  increases  at  the  successive  levels  through  the  GS  structure.  Thus 
the  cancellation  process  in  the  main  channel  is  terminated  at  the  level  where  the  noise  power  is  mea¬ 
sured  to  be  a  minimum.  Hence  the  number  of  levels  prior  to  this  termination  point  is  estimated  to  be 
equal  to  ND oF.  If  we  compute  AT '3d B  for  the  DOF’s  monitoring  technique,  then 


^3dB  -  2Ndof  -  2. 


As  a  result  for  the  three  techniques 

normal  : 

noise 
injection  ‘ 

DOF 

monitoring  ‘ 


N 

iJaux 


*3dB 

N 

Iwawc 


N 

l'aux 


=  2 


«  1  + 


N  DOF 
^ aux 


-  WpQF 
N 

aux 


(90) 

(91) 


(92) 


(93) 


where  Naux  =  N  -  1 . 

We  plot  these  quantities  vs  Ndof^u  in  Fig.  7.  Note  that  the  DOF  monitoring  technique  is 
faster  converging  than  is  the  noise  injection  technique.  Hence,  DOF  monitoring  with  the  GS  can¬ 
celler  implementation  offers  a  computationally  efficient  and  fast  converging  adaptive  canceller  algo¬ 
rithm.  Also,  because  B  is  a  measure  of  how  perturbed  the  adaptive  weights  are  about  their 
optimal  values,  we  see  that  the  transient  sidelobe  settling  times  should  be  faster  by  the  use  of  DOF 
monitoring. 

IX.  OTHER  SIDELOBE  LEVEL  REDUCTION  TECHNIQUES 

In  this  section,  we  briefly  discuss  other  techniques  that  can  be  used  to  reduce  the  average  tran¬ 
sient  sidelobe  level.  The  basis  of  these  techniques  is  to  reduce  the  number  of  auxiliar /  input  channels 
N  -  1  of  the  SLC  to  the  necessary  number  of  DOFs,  N^qf- 
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Fig.  7  —  Kmb/N—x  vs  Noor/Nmu  for  three  SLC  techniques 

The  first  technique  is  to  form  N  -  1  orthogonal  auxiliary  beams  that  cover  the  0  to  2x  angular 
space.  This  can  be  done  by  setting  the  unitary  matrix  transform  U  equal  to  the  Butler  matrix  P 
where 


and 


v.  j 


n,  m  =  1,2,  •  •• ,  N  -  1 


(94) 


—  •  2y 

r*_,  =  e  ,  j  =  V^T.  (95) 

Thereafter,  only  those  beams  that  contain  high  levels  of  external  interference  are  inputted  into  the 
SLC  SMI  canceller.  Hence  the  number  of  nonsignificant  eigenvalues  is  reduced  so  that  A SLa  will  be 
reduced.  Also,  the  effects  on  the  main  channel’s  sidelobe  level  are  localized  around  where  (in  angle) 
the  various  external  interferences  appear  .  Problems  involved  with  this  technique  are  the  significant 
amount  of  hardware  necessary  to  generate  N  -  1  beams  and  the  assurance  that  enough  beams  are 
used  to  suppress  the  externa,  noise  environment  (there  may  be  multiple  interferences  in  a  single 
beam). 

A  second  technique  related  to  the  first  technique  is  to  generate  auxiliary  beams  that  point  at  the 
various  external  interferences.  Even  though  the  auxiliary  beams  may  not  be  orthogonal,  the  effects 
on  the  main  channel’s  sidelobe  level  are  localized  around  where  in  angle  the  various  interferences 
appear.  A  prerequisite  for  this  technique  is  that  the  interference’s  approximate  (within  an  auxiliary 
antenna  beamwidth)  angular  position  be  known.  This  can  be  done  by  use  of  a  variety  of  techniques 
(for  example,  superresolution  [10]). 
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Appendix 

UNITARY  MATRIX  TRANSFORM  OF  A  STEERING  VECTOR 

In  this  appendix,  we  give  a  simple  procedure  for  generating  a  unitary  matrix  transform  U, 
which  transforms  the  following  steering  vector 

s  *  ^(1.  . 'W-mf  (A1) 

into  the  (1,  0,  0,  ... ,  0)r  vector.  First,  we  define  an  N  x  N  diagonal  matrix  5,  such  that  the  diag¬ 
onal  elements  are  equal  to  the  corresponding  elements  of  the  steering  vector  s.  Next,  we  select  any 
unitary  transformation  U0,  such  that  all  of  the  first  row  elements  equal  1/Vn.  Examples  of  this  kind 
of  matrix  are  the  Butler  (see  Eq.  (90))  and  the  Hadamard  matrices.  Finally  we  set 

U  «  UqS*  .  (A2) 

It  is  elementary  to  verify  that  UU'  =  /  and  that  UqS*s  =  (1,  0,  0 . 0)T. 
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